** Replication files for "Fuel on the Fire: Has the Politics of COVID-19 Accelerated Secularization in America?"


** Note: you will need to replace "ADD FOLDER" with the correct information for your files **

** Study 2 (first experiment) **

** use "<ADD FOLDER> Fuel on the Fire, Experiment 1.dta"

set level 95

** Figure 1; Table A4 **
regress personsec_add_01_nosci3 covid_rel2 covid_norel2 if pid_3==1 & correct_answer2>0
estimates store sec_dems


margins, asobserved at(covid_rel2=0 covid_norel2=0)
margins, asobserved at(covid_rel2=1 covid_norel2=0)
margins, asobserved at(covid_rel2=0 covid_norel2=1)


regress personsec_add_01_nosci3 covid_rel2 covid_norel2 if pid_3==3 & correct_answer2 >0
estimates store sec_reps


margins, asobserved at(covid_rel2=0 covid_norel2=0)
margins, asobserved at(covid_rel2=1 covid_norel2=0)
margins, asobserved at(covid_rel2=0 covid_norel2=1)

etable, estimates(sec_dems sec_reps) mstat(N) mstat(r2_a) showstars showstarsnote export("<ADD FOLDER>", replace)


coefplot (sec_dems, label(Democrats)) (sec_reps, label(Republicans) msymbol(T)), drop(_cons) xline(0) levels(95) title("Effect on Secularism")
graph save  "<ADD FOLDER>", replace


** Table A5 **
** without attention check **

regress personsec_add_01_nosci3 covid_rel2 covid_norel2 if pid_3==1 
estimates store sec_dems_nocheck

regress personsec_add_01_nosci3 covid_rel2 covid_norel2 if pid_3==3
estimates store sec_reps_nocheck

etable, estimates(sec_dems_nocheck sec_reps_nocheck) mstat(N) mstat(r2_a) showstars showstarsnote export("<ADD FOLDER>", replace)


** Study 3 (second experiment) **

** use "<ADD FOLDER"> Ful on the Fire, Experiment 2.dta"

**Table A6 **

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==1 & correct_answer2 >0
estimates store sec_dems_exp2

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==3 & correct_answer2 >0
estimates store sec_reps_exp2


etable, estimates(sec_dems_exp2 sec_reps_exp2) mstat(N) mstat(r2_a) showstars showstarsnote export("<ADD FOLDER>", replace)


** Table A7 **
** without attention check **

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==1 
estimates store sec_dems_exp2_nocheck

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==3 
estimates store sec_reps_exp2_nocheck


etable, estimates(sec_dems_exp2_nocheck sec_reps_exp2_nocheck) mstat(N) mstat(r2_a) showstars showstarsnote export("<ADD FOLDER>", replace)


** Figure 2, upper panel; Table A8 ** 

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==1 & vaccines <4 & correct_answer2 >0
estimates store sec_dems_novax_exp2


margins, asobserved at(covid_rel=0 covid_norel=0)
margins, asobserved at(covid_rel=1 covid_norel=0)
margins, asobserved at(covid_rel=0 covid_norel=1)

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==1 & vaccines >3 & correct_answer2 >0
estimates store sec_dems_provax_exp2

margins, asobserved at(covid_rel=0 covid_norel=0)
margins, asobserved at(covid_rel=1 covid_norel=0)
margins, asobserved at(covid_rel=0 covid_norel=1)


etable, estimates(sec_dems_novax_exp2 sec_dems_provax_exp2 ) mstat(N) mstat(r2_a) showstars showstarsnote export("<ADD FOLDER>", replace)


coefplot (sec_dems_provax_exp2, label(Pro-Vaccine Democrats)) (sec_dems_novax_exp2, label(Anti-Vaccine Democrats) msymbol(T)), drop(_cons) xline(0) levels(95) title("Effect on Secularism")
graph save "<ADD FOLDER>", replace


** Table A9 **
** without attention check **

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==1 & vaccines <4 
estimates store sec_dems_novax_2_nocheck

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==1 & vaccines >3
estimates store sec_dems_provax_2_nocheck

etable, estimates(sec_dems_novax_2_nocheck sec_dems_provax_2_nocheck ) mstat(N) mstat(r2_a) showstars showstarsnote export("/Users/dcampbe4/Dropbox/Secularism Project/Secular Surge II papers/covid paper/PB revision/table_a3_nocheck.docx", replace)


** Figure 2, lower panel; Table A10 **

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==3 & vaccines <4 & correct_answer2 >0
estimates store sec_reps_2_novax

margins, asobserved at(covid_rel=0 covid_norel=0)
margins, asobserved at(covid_rel=1 covid_norel=0)
margins, asobserved at(covid_rel=0 covid_norel=1)

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==3 & vaccines >3 & correct_answer2 >0
estimates store sec_reps_2_provax

margins, asobserved at(covid_rel=0 covid_norel=0)
margins, asobserved at(covid_rel=1 covid_norel=0)
margins, asobserved at(covid_rel=0 covid_norel=1)


etable, estimates(sec_reps_2_novax sec_reps_2_provax) mstat(N) mstat(r2_a) showstars showstarsnote export("<ADD FOLDER>", replace)


coefplot (sec_reps_2_provax, label(Pro-Vaccine Republicans)) (sec_reps_2_novax, label(Anti-Vaccine Republicans) msymbol(T)), drop(_cons) xline(0) levels(95) title("Effect on Secularism")
graph save  "<ADD FOLDER", replace


** Table A11 **
** without attention check **

regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==3 & vaccines <4 
estimates store sec_reps_2_novax_nocheck


regress personsec_add_01_nosci3  covid_rel covid_norel if pid_3==3 & vaccines >3 
estimates store sec_reps_2_provax_nocheck

etable, estimates(sec_reps_2_novax_nocheck sec_reps_2_provax_nocheck) mstat(N) mstat(r2_a) showstars showstarsnote export("<ADD FOLDER>", replace)

